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Abstract. We consider integer programming and the semi-group 
membership problem. We develop and extend the approach started 
in [5] [6] so as to provide the following theorem of the alternative: 
the system b = Ax has no nonnegative integral solution x E N n if 
and only if p(b) < for some given polynomial p. The coefficients 
of p form a vector which lies in some convex cone Q, and so we 
characterize f2. We also provide a hierarchy of linear programming 
relaxations, where the continuous case Ax = b describes the first 
relaxation in the hierarchy for x E M. n and x > 0. 



1. Introduction 

This paper is concerned with certificates for integer programming 
(IP) as well as with the semi-group membership problem. That is, 
given a finitely generated abelian group G (e.g. Z m ), a semi-group 
G a C G generated by a finite family (afe)^ =1 C G, and an element 
b G G, we provide a certificate of b G G a or b G" G a . We build upon and 
extend previous work of [51 El E], notably on a discrete Farkas lemma 
for IP. Among other things, we provide a hierarchy of linear program- 
ming relaxations (LP-relaxations) for integer programming. The first 
relaxation in the hierarchy is just the usual LP relaxation, which then 
appears as a first-order (or linear) approximation to the discrete case, 
whereas usually IP is viewed as an arithmetic refinement of LP. We 
also provide a theorem of the alternative (or duality theorem) in the 
form of a polynomial certificate associated with the IP problem, and 
we compare with the certificate for LP obtained by the standard Farkas 
lemma. 

A central idea in nonconvex optimization is to replace a non convex 
(hence hard) problem with a suitable easier convex problem in some 
lifted space, but at the price of increasing the dimension. For instance 
in the lift-and-project approach for polynomial optimization (e.g. 0- 
1 problems) one replaces x G W 1 with the vector y = (x a ) of all 
moments and solves some hierarchy of appropriate linear or semidefinite 
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relaxations. The interested reader is referred for more details to e.g. 
Sherali and Adams [121 H2], Lovasz and Schrijver [TTJ, and Lasserre 
[7J[8]; see also Laurent [10] for a comparison. Of course, any IP problem 
can also be modeled via polynomial equations. For example, if the entry 
Xi is bounded by some integer M, then one may include the polynomial 
constraint nfclo( x * ~~ ^) which forces Xj to be an integer, and so the 
above methodology applies. However, since the degree in the constraint 
is M (as opposed to 2 in the Boolean case), the size of the first linear 
or semidefinite relaxation in the hierarchy (in e.g. [Sj) is already very 
large because it includes moments up to order M. 

Let A G pf" xn De a fixed matrix. In the approach developed in [51 [6] 
the integral solution i 6 N" to Ax = b is also lifted to some y G MP 
with p < n Ylj(l+bj). But this time the lifting process has a different 
meaning. Indeed there is some very simple matrix E G R nxp such that 
x := Ey G R™ is now a point in the integer hull of feasible solutions. 
Furthermore, the vector y is a point of some polytope and several 
interpretations can be deduced from the lifting process. For example it 
was already shown in [5] that the lifting process can be used to prove 
that b = Ax for some nonnegative integral vector x G N™ if and only if 
the polynomial z b — 1 has a nonnegative representation in the binomial 
ideal generated by the binomials (z Ak — 1). But this interpretation is 
only one among many other interpretations as shown in the present 
paper. 

Interestingly, one can also use the lifting process to provide a hierar- 
chy of LP relaxations for the IP problem, so that the continuous case 
{Ax = 6:xGlR n ,x>0} appears as a first-order (or linear) approxi- 
mation of the discrete case with x G N n . This hierarchy also provides 
us with a theorem of the alternative (or duality theorem) which uses 
a non linear polynomial, a discrete analogue of the celebrated Farkas 
Lemma in linear algebra and convex optimization. Recall that the 
Farkas Lemma provides a membership certificate for the convex cone 
O := {Ax : x G lR n , x > 0}, in the form b (jL if and only if u'b < for 
some uj with u'A > 0. In other words b ^ O if and only if the linear 
polynomial z t—> p u (z) := uj'z is negative when evaluated at z — b. 



Contribution. Every finitely generated abelian group G can be iden- 
tified with a subset of Z m and so: 

• We firstly show that the semi-group membership problem reduces 
to the existence of a nonnegative integral vector x G N n , solution of 
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some related linear system Ax = b for some nonnegative integral ma- 
trix A G N mxn . Shevshenko [TBI, p. 11] developed a similar result in 
the framework of additive semigroups with unity, contained in a finite 
generated abelian group. 

• Set p := Yuk=i TYjLii^+bj— Aj. k ). We next show that existence of 
a nonnegative integral solution x G N n to the linear system Ax = b 
reduces to the existence of a nonnegative real solution y G MP for a 
system of linear equations of the form : 

(1.1) bj = y'A^\ Vl<j<m; 

(1.2) bibj = y 'A iei+e i\ V 1 < i < j < to; 



VO < Zj < bj, 
z x A Yz m = 5; 



(1.3) &T = y'A^, 

(1.4) #---«fr = y'A®; 

for some appropriate nonnegative integer vectors A^ ej \ A^ ei+Sj \ A^ G 
W with z G N m and z < b. The parameter 5 > 1 in (11.31) is the 
degree of the monomial b i— > fe^ 1 ■ ■ • 6^ in (II. 3p . Therefore a certificate 
of b ytz Ax for every x G N n is obtained as soon as any subsystem 
of (ll.ip - (ll.4l) has no solution y G M p ; that is, one does not need to 
consider the entire system (ll.ip - (ll.4l) . 

We can index the entries of A^ = (A^[k, u}) G W in (11.11) in such 
a way that A^ e ^ [k, u] = Aj.k for all u, 1 < j < m, and 1 < k < n. If we 
also index the entries of y = (y[k, u}) G M? in the same way, the new 
vector x = (xk) G lR n with Xk ■= J2 U y[k, u] satisfies Ax = b and belongs 
to the integer hull of {i 6 I" : Ax = b, x > 0}, whenever y is a solution 
to (ll.ll) - (ll.4p . In this approach, LP (or the continuous case) appears 
as a particular "first order" (or "linear") approximation of IP (the 
discrete case). Indeed, if one considers (11.11) alone (i.e. ignoring (11.21) - 
(11.41) which have nonlinear right-hand-sides terms b z for ||z|| > 1) then 
from any nonnegative solution y of (II. ip one obtains a real nonnegative 
solution x G M n of Ax = b, and conversely. 

To construct a natural hierarchy of LP-relaxations for the IP feasi- 
bility problem Ax = b, x G N n , just consider an increasing number of 
equations among the system (ll.ll) - (ll.4p . so that the last (and largest 
size) LP-relaxation is the whole system (ll.ip - (ll.4l) that describes the 
integer hull of the set {x G N n : Ax = b}. Thus, if on the one hand the 
discrete case is an arithmetic refinement of the continuous one, on the 
other hand the discrete case can be approximated via LP-relaxations of 
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increasing sizes, and these relaxations are different from the lift-and- 
project ones described in e.g. ([ID])- To the best of our knowledge 
such a hierarchy has not been investigated before. Even if it is not 
clear at the moment whether this hierarchy of linear relaxations is use- 
ful from a computational viewpoint, it provides new insights for integer 
programming. 

On the other hand it was already proved in j5j [6] that existence of 
a nonnegative integral solution x G W 1 to the linear system Ax = b, 
reduces to the existence of a nonnegative real solution y G M? for a 
system of the form 

(1.5) (-1,0,. ..,0,1)' = 0y, 

where is some appropriated network matrix (hence totally unimodu- 
lar). We show that the system (|l.l[) - (|1.4p can be deduced from (jl.5p by 
multiplying it from the left times a square invertible matrix A G ~R SXS . 
In particular A is a Kronecker product of Vandermonde matrices. 
Hence any real vector y > is solution of (11.51) if and only the same y 
is solution of (OI)- (Oj) . 

• We provide a polyhedral convex cone Oct 5 associated with (11.11) - 
(11.41) for some s6N, such that a direct application of Farkas lemma to 
the continuous system (ll.ip - (ll.4l) implies that either b = Ax for some 
integral vector x G N n or there exists £ = (£„,) G Q such that p^(b) < 
for a polynomial p^ G M.[ui, u m ) of the form 

(1.6) u^p^u) = ^ uW - 

Thus (11.61) provides an explicit nonlinear polynomial certificate for IP, 
in contrast with the linear polynomial certificate fro LP obtained from 
the classical Farkas lemma. 

In the discrete Farkas lemma presented in [51 E] the author defines a 
polyhedral cone Q 2 C M s associated with (11.51) . and proves that either 
b = Ax for some nonnegative integral vector x G N n or there exists 
7i G O2 such that (—1, 0, 0, 1) ■ it < 0. It turns out that the vector 
£ = (£«>) i n (11-61) indeed satisfies 71 = A'£ for a square invertible matrix 
A defined as the Kronecker product of Vandermonde matrices. 

• Inspired by the relationships between the systems (ll.ip - (ll.4l) and 
(11.51) . we finally show that existence of a nonnegative integral solution 
x G N n for the linear system Ax = b reduces to the existence of a 
nonnegative real solution y G MP for a system of linear equations of the 
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form: 

M (-1,0, ...,0,1)' = M6y, 

where M G IR SXS is any square invertible matrix. We can apply the 
standard Farkas lemma to any one of the linear systems presented 
above, and deduce a specific (Farkas) certificate for each choice the 
invertible matrix M. Each certificate can be seen as a theorem of the 
alternative of the form: the system b = Ax has no nonnegative integral 
solution x G N n if and only if f(b) < for some given function /. 
For the particular choice M := A', existence of a nonnegative integral 
solution x G N n for the linear system Ax = b reduces to existence of 
a nonnegative real solution y G MP for a system of linear equations of 
the form: 

u b - 1 = y'D [u \ V u G N m with u < b, 

for some appropriate nonnegative integer vectors G N p . In this 
case, the function / involved in the Farkas certificate f(b) < has the 
exponential-like expansion 

(1.7) f(u) := (*"-!)■ 

zeN m 

Both certificates (11.61) and (jl.7p are different from the certificate 
obtained from the superadditive dual approach of Gomory and Johnson 
[2], Johnson [4], and Wolsey [IB] . In particular, in [TBJ the linear system 
associated with such a certificate has dimension s 2 x s, which is larger 
than the size ns x s of the system (11.51) of this paper. 

The method. A theorem of the alternative is obtained in three steps: 

{1} One first shows that a linear system Ax = b (with A G N mxn 
and b G N m ) has a nonnegative integral solution x G N n if and only 
if the function / : Z m — >M. given by z i— > /(z) := 6 2 can be written as 
a linear combination of functions z i— > := (m+A^) 2 — m z weighted 

by some nonnegative coefficients (for the indexes u G N m with u <b). 

{2} One then shows that computing the nonnegative coefficients in 
the above linear decomposition of / is equivalent to finding a nonnega- 
tive real solution to a finite system of linear equations whose dimension 
is bounded from above by ns x s with s := rij(l+^')- 

{3} Finally one applies the standard continuous Farkas lemma to 
the linear system described in {2} and obtains the certificate (11.61) . 

This approach is similar in flavor but different from the one in [51 [6J . 
However, they are strictly equivalent and can be deduced one from each 
other under some appropriate linear transformation whose associated 
matrix A G N sxs has a simple explicit form. 
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2. Notation and definitions 

The notation R, Z and N = {0, 1,2, ...} stand for the usual sets of 
real, integer and natural numbers, respectively. Moreover, the set of 
positive integer numbers is denoted by N* = N\ {0}. Given any vector 
b G Z m and matrix A G Z mxn , the [/c]-entry of b (resp. \j; fc]-element of 
A) is denoted by either 6 fc or b[k] (resp. Aj. k or A\j\ k]). The notation 
A' stands for the transpose of any matrix (or vector) A G M mxn ; and 
the fcth column of the matrix A is denoted by A k := (v4 1;fc , A m . k )' . 



The semi-group membership problem. A classical result in group 
theory states that every abelian group G with m G N* generators is 
isomorphic to some Cartesian product 



m—ql 



G = [Z/piZ] x [Z/p 2 Z] x ■ ■ ■ x [Z/p,Z] x [Z 

for some set of numbers {pj} C N* and q < m; see e.g. pQ. One may 
even suppose that every pj divides p k whenever j < k. Therefore, if 
one introduces the extended m-dimensional vector: 

(2.1) P := (pi,p2,...,p„oo,...,oo)', 

the abelian group G is isomorphic to the group G of vectors x e Z m 
such that < Xj < pj for every 1 < j < q; notice that q < m. The 
group sum x © y of two elements x and j/inGc Z m is then defined by 

(2.2) x © y := (x + y) mod P in G, 

where the sum x + y is the standard addition on Z m and the modulus 
of the sum (x + y) mod P is calculated entry by entry, so that for every 
index 1 < k < m, 



(2.3) 



[(x + y) modP] fc = | 



{xk + Vk) mod P k if Pfc < oo, 
(x k + yk) if Pk = oo. 



Hence from now on we suppose that G = G C Z m . Next let {a k } C G 
be a collection of n elements of G. Each element ak can be seen as a 
vector of Z m , for 1 < k < n, so that the semi-group generated by {ak} 
is the same as the set 

G a := {Ax mod P\x <E N m }, with 
A := [ ai \a 2 \---\a n }eZ mxn . 

Thus, given b G G, the semi-group membership problem of deciding 
whether b G G a is equivalent to deciding whether the system of linear 
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equations b = Ax mod P has a solution x G N n . This in turn is 
equivalent to deciding whether the following system of linear equations 

-B B\ (u 
I ' [w. 



(2.4) b = Ax + 



has a solution (x,u,w) in N™ x f x N q , where 
(2.5) B : 



(p x ••• 0\ 
p 2 ■■■ 



g N 9X9 . 



\o ■■■ pj 

Hence, with no loss of generality, the membership problem is equiva- 
lent to deciding whether some related system of linear equations Ax = b 
(with A G Z mx£ and b G Z m ) has a solution x G N £ , which is the prob- 
lem we will consider in the sequel. 

3. Existence of integer solutions to {Ax = b, x > 0} 

Let WL[z] = M[zi, . . . , z m ] be the ring of real polynomials in the vari- 
ables z = (z\, . . . , z m ) of W 71 . With A G Z mxi we analyze existence of 
a nonnegative integer solution x G N e to the system of linear equations 
Ax = b. 

According to the computational complexity terminology in [14J, let 
ip be the facet complexity of the rational convex polyhedron P := {x G 
R e : x > 0,Ax = b}. That is, each inequality in Ax < b, Ax > b, 
and x > has size at most (p. Corollary 17.1b in [TH p. 239] states 
that P contains an integral vector of size at most 6£ 3 (p, if P contains 
an integral vector x G W. Hence existence of an integral vector in P 
is equivalent to analyze existence of an integral vector of the convex 
(compact) polytope 



{x G R e : x > 0, Ax = b, < M A>b }, 



i=i 



where Ma^ is obtained explicitly from the facet complexity of P. 



3.1. Reduction to A G N mxn . In view of the above one may restrict 
our analysis to the existence of a nonnegative integral solution x G W 1 
for a system of linear equations Ax = b associated with a rational 
convex polytope P := {i 6 K" : Ax = b, x > 0}, where A G Z mxn and 
b G Z m . But this in turn implies that one may restrict our analysis to 
existence of a nonnegative integral solution y G N n+1 to a system of 
linear equations A*y = b* where A* G N[ m + 1 ] x [™+ 1 ] and b* G N m+1 . 
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Indeed, if A G Z mxn and A £ N mxn , let a G N n be such that 

(3.1) Aj. k := Aj. k + a k >0; V 1 < j < m, 1 < k < n. 
Since P is a (compact) polytope, 

(3.2) p := max < > a k x k > < oo. 

k k=i J 

In particular p G N. Now define b G N m by 

(3.3) b:=b + pe m >0 with e m := (1, 1)' G N m . 

Let A G N mxn be defined as in (jOl . The solutions x G N n to 
the original system Ax = 6 are in one-to-one correspondence with the 
solutions (x, u) G W 1 x N to the extended system 

/o ^\ Ax + e m M = 6, 

a'x + tt = p. 

Indeed, if Ax = b with x G N n , then 



Ax + e r 



n 



k=l 



- e m (a'x) + pe m = b + pe r 



The following identity follows from the definitions for A G N mxn and 
b G N m given in (j3TT]) -([3T3l : 

Ax + e m u = b with u := p — a'x G Z. 

Notice that u > because p > a'x according to (13.21) . so that (x,u) 
is the integer nonnegative solution of (13 .4p that we are looking for. 
Conversely let (x,u) G Z n+1 be a solution to (13.41) . The definitions for 
A G N mxn and b G N m given in (pnj) - (pO]) imply that 



jA.X ~\~ Qf> 



n 



E 



ttfcXfc 



+ e m w = 6 + p e r , 



k=l 

so that Ax = b with x G N n because u = p—a'x. Hence the existence 
any solution x G N n to Ax = b is completely equivalent to the existence 
any solution (x, u) G N n+1 to 



b ,r ■ r i with a* ^ 4 



7)1 



p J \u J \a' 1 

and the new matrix A* G ^[ m + 1 ] x [ n + 1 ] has only nonnegative integer 
entries. 
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3.2. The main result. Given vectors b £ N m and z £ M. m , the nota- 
tion z b stands for the monomial z^z^ ' ' ' z rn We also need to 
define a pair of matrices A and that we will use in the sequel. 

Definition 1. Let A £ N mxn and (3 £ N m be such that A k < (3 for 
each index 1 < k < n. Set the integers 

m n m 

(3.5) s : = n( 1+ ^') and P=EII i 1 +^- A r,k) < ns. 

j=l k=l j=l 

Let A £ N sxs be a square Vandermonde matrix whose rows and 
columns are indexed with the nonnegative vectors z,w £ N m (according 
to e.g. the lexicographic ordering) so that the [z;w]- entry is given by 

(3.6) A[z;w] =w z = w z 1 1 w z 2 2 ---w z ™ V z,w < (3; 

and where we use the convention that 0° = 1. 

Let £ Z sxp be a network matrix whose rows and columns are 
respectively indexed with the nonnegative vectors w £ N m and (k, u) £ 
NxN m ; so that the (w; (k,u))- entry is given by 

{-1 ifw = u<(3-A k , 
1 if w = U +A k < (3, 
otherwise, 

for all indexes w < (3, u < (3—A k , and 1 < k < n. 

The following result is straightforward. 

Lemma 2. The square Vandermonde matrix A £ W xs in Ii3. 6\) is 
invertible. The matrix £ ~R sxp in ^3. 7| ) is a network matrix, and so 
it is totally unimodular. 

Proof. It is easy to see that A is invertible because it is the Kroneker 
product £>W ® <g> • • • ® IlH of m square Vandermonde matrices 



(3.8) D [j] 



fl 1° 2° ••• {(3 3 f\ 
I 1 2 1 ... (/3,) 1 



\0 1^ 2^ ••■ 

see e.g. (21 HZj. And each is obviously invertible. By inspection 
it turns out that is a network matrix, that is, it is a matrix with 
only {0, ±1} entries and with exactly two nonzero entries 1 and —1 in 
each column. Therefore is totally unimodular; see e.g. Schrijver [TU 
p. 274]. □ 
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Definition 3. Given a finite set U C N m and a collection of coefficients 
q[u] G R ; we say that the polynomial 

has multivariate degree bounded by a vector [3 G N m if and only ifu<(3 
for every u G U . 

The following technical result was essentially shown in (5j E] but we 
include the proof for the sake of completeness. 

Theorem 4. Let (3,beW l and A G N mxn be such that : 

(3.9) (3 > b and (3 > A k for all 1 < k < n. 

The following three statements (a), (b) and (c) are all equivalent : 

(a) : The linear system Ax = b has a solution x G N n . 

(b) : The polynomial z i— > z b — 1 := z^z^ 2 ■ • • z h ™ — \ can be written 
as follows : 

n 

(3.10) z b -l = Y,Qk(z)(z Ak -l) 

k=l 

for some real polynomials Qk G M.[zi, z 2 , z m ] with nonnega- 
tive coefficients and multivariate degree bounded by the vector 
f3—Ak, for all 1 < k < n. 

(c) : There is a real nonnegative solution y G MP for the system of 
linear equations 

(3.11) b = 6y, 

where G Z sxp is given as in 7p of Definition^ and the new 
vector b G If has entries indexed with the nonnegative vector 
w G N m , so that the [w]-entry is given by 

( -1 if w = 0, 

(3.12) bH := < 1 ifw = b, V0<w<(3. 

[ otherwise. 

Notice that 6 Z sxp in (13.71) and (13. lip is determined only by the 
entries of (3 G N m and A G N mx ", so that 6 is independent of b G 7L m . 
Observe that if b > A^ for every 1 < k < n, then one may take (3 := b. 
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Proof, (a) (b). Suppose that b = Ax for some x G N n . Consider 



XI— 1 

E ^ 

g=0 



the following polynomials (where we use the convention (-) = 0) 

XI — 1 

z h-> Qi(z) 
z i-> Q 2 (z) 
z i-> Q 3 (z) 



nA x 



q=0 



ra-1 

n 

A;=i 



X3-I 

E* 

<?=0 
n In- 1 



9^3 



,XkAk 



E 2 

9 =0 



It is easy to see that the polynomials Qk satisfy equation (I3.10p . 
Moreover each Qk has nonnegative coefficients, and the multivariate 
degree of Qk is bounded by the vector (3— Ak for 1 < k < n because 
b<(3. 

(b) <^> (c). Existence of polynomials Qk £ R.[z\, Z2, z m ] with 
nonnegative coefficients and multivariate degree bounded by the vec- 
tor (3—Ak (for all 1 < k < n) is equivalent to the existence of real 
coefficients u) > such that 

Qk(z)= v[k,u}z u v l<A;<n. 

u& m ,u<l3-A k 

Then rewrite equation (13 . 1 j) as follows: 



(3.13) 



1 = E E y^ u K zU 



k=l ueN m , 
u<fi-A h 



The vector of coefficients y = (y[k,u\) > satisfies (13.131) if and only 
if it satisfies the following system of linear equations 



(3.14) 



/y[l,(0,...,0,0)]\ 
(0, .... 0, 1)] 

\y[n, (P-A n )]J 
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where b 6 Z s is given in ( I3.12p and the matrix G ^ sx p j s given in 
(13.71) . that is, such that 

{—1 if w = u < b—Ak, 
1 if w = u+A k < b, 
otherwise. 

Each row in (13.141) is indexed with the monomial z w G N m , for w G 
N m and w < b (according to e.g. the lexicographic ordering). 

(c) =>• (a). Let y = (y[k, u)) > be a real vector such that (13.111) and 
( 13. 14[) hold. The matrix is totally unimodular according to Lemma [2} 
and so there exists a nonnegative integer solution y = (y[k, u\) > to 
( 13. lip and ( 13.141) because the left-hand-side of ( 13.141) is an integral 
vector. Therefore (I3.13P also holds with the new vector y, that is: 

n 

k=i «eN m , 

u<P-A k 

Given a fixed index 1 < j ' < m, differentiate both sides of the equation 
above with respect to the variable Zj and evaluate at the point z = 
(1, 1, 1), so as to obtain: 

n 

(3.15) bj = ^Aj.kXk, with x k := ^2 v[hu). 

k=l ueN m , 

The nonnegative integer vector x = (xi, ...,x n )' G N m satisfies the 
desired result Ax = b. □ 

Remark 5. One linear constraint in is redundant, because the 

addition of all the rows in A3.11\) - f3.1fy) yields the trivial equality = 
Oy (recall that is a matrix with only {0, ±1} entries and with exactly 
two nonzero entries 1 and —1 in each column). 

Remark 6. The matrix G Z sxp in ^3. 7\ ) is independent of b and 
contains all information about all the integer problems Ax = b for 
x G N™ and be N m with b<(3. 

The linear system b = y in (13. lip of Theorem H] is quite interesting 
in the sense that existence of a real solution y > is totally equiva- 
lent to the existence of a nonnegative integer solution igN" for the 
problem Ax = b, but it is not easy to see at first glance what is the 
relation between the two solutions y and x. Moreover equation (13.111) 
in Theorem H] is not unique at all. We can multiply both sides of (13.111) 
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by any square invertible matrix M G W xs and we obtain that y > is 
a real solution to (13. Ill) if and only if y is also a solution for 

Mb = M9y. 

An interesting issue is to determine what is an adequate matrix M 
such that relationship between the solutions y > and x G N n is 
evident. We claim that one of such matrix M is the square Vander- 
monde matrix A G N sxs presented in (I3.6P of Definition [TJ We use the 
convention that 0° = 1 and 2 = whenever ^ z G N m . 

Theorem 7. Let (3,beN m and A G N mxn be such that : 

(3.16) (3>b and (3 > A k for all 1 < k < n. 

The following three statements (a), (b) and (c) are all equivalent : 

(a) : The linear system Ax = b has a solution igN". 

(b) : There is a real nonnegative solution y = (y[k,u\) G W to 
the system of linear equations 

(3.17) b z -o z = J2 E y[k,u]{(u+A k y-u z ), 

k=l «eN m , 
u</3-A k 

for every z G N m with < z < (3. 

(c) : There is a real nonnegative solution y = (y[k, u]) G W to the 
system of linear equations 

n 

(3.18) ^-1 = E E y[k,u}{z u+A *-z u ), 

k=l ueN m , 

for every z G N m with < z < (3. 

Proof. The equivalence (a) <^ (b) easily follows from Theorem S] after 
noticing that the [z; (k, w)]-entry (resp. [^]-entry) of the product A© 
(resp. Ab) is given by 

(3.19) (AQ)[z;(k,u)) = ((u+A k ) z — u z ) G N, 

(3.20) resp. (Ab)[z\ = (b z -0 z ), 

for all nonnegative indexes z G N m and (k, u) G MxN m with z < (3, 
u < (3—A k and 1 < k < n. 
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The equivalence (a) <^ (c) is proved in a similar way; we only need 
to observe that multiplying by the transpose A' is equivalent to inter- 
change the exponentials z — > w z by powers z — > z w , so that: 

(A'e)[z;(k,u)} = (z u+Ak - z 11 ) eN, 
and (A'b)[z] = (z b -z°). 

□ 

Remark 8. A nonnegative real vector y eMP is a solution of equation 
113.11]) in Theorem if and only if y is also a solution of $3.17\ ) and 



$3.18]) . Moreover the linear system $3.18]) can be directly deduced from 
$3.10]) by evaluating at z G N m with z < b. 



4. A HIERARCHY OF LINEAR PROGRAMMING RELAXATIONS 

Let 0, b G N m and A G N mxn be such that (3 > b and (3 > A k for all 
1 < k < n. Consider the matrices A G M sxs and G M sxp given in 
equations (13.6p - (13.7p of Definition [TJ and define the new vector b G Z s 
whose entries are indexed with the non negative vector w G N m , so 
that the [fc]-entry is given by (13.121) : 

( -1 ifw = 0, 

(4.1) bH := < 1 if iu = 6, V0<w<(3. 

[ otherwise, 

The fact that A is invertible (according to Lemma [2]) allows us to 
give two equivalent definitions of the following polytope in MP, 

(4.2) V := {y G W : y > 0, A6y = Ab} 

= {y G W : y > 0, 0y = b} C MP. 

As it is state in Remark [HJ any nonnegative vector y G MP is a 
solution to (13.171) if and only if y lies in V . Let y = (y[k, u]) be an 
element of V. The entries of y are indexed according to Theorems H] 
and [71 Equation (13.171) implies that x G lR n with x k := ^2 u y[k,u], is 
a solution to Ax = b. Indeed bj = J2k ^j;k x k for 1 < j < m after 
evaluating (I3.17P at z = ej, the basic vector whose entries are all equal 
to zero, except the j-entry which is equal to one, so that w ej = Wj for all 
w G Z m . The relationship between x and y comes from a multiplication 
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by some matrix E G N nxp , that is : 



(4.3) 



x 



Ey with E :- 




\0...0 0...0 




1...1J 



where p k := Yl™ =1 (l+/3j—Aj. k ) for 1 < k < n, so that p = J2 k Pk 
according to (13. 5p in Definition [IJ We obtain the following result as a 
consequence of Theorem [7J 

Corollary 9. Let (3, b G N m and A e N mxn 6e snc/i tfmt (3 > b and 
(3 > Ah for each index 1 < k < n. 

(a): Existence of a nonnegative integral solution x G W 1 to the linear 
system Ax = b reduces to the existence of a nonnegative real solution 
y = (y[k,u\) G M? to the system of linear equations : 

(4.4) b t = y'A^\ V 1 < j < m; 

(4.5) bibj = y 'A {ei+ej \ V 1 < % < j < m; 



(4.6) 



b zi ■ ■ ■ h Zm 



r'A {z \ 



VO < zj < b d , 
Z\-\- ■ ■ ■ -\-z m = 5: 



(4.7) b^---bt = y'A^- 

for some appropriate nonnegative integer vectors A^\ A^ £l+ej \ A^ G 
N p with z G N m and z < (3. The parameter 8 > 1 in l[4-6\ ) is the degree 
of the monomial b i— ► b z , and we use the indexation for the entries 
°f y given in Theorems [7] and [?| In particular all vectors A^ are 
independent ofb, and the entries of A^ e ^ G W in ^4-4\) are given by : 



A^[k, 



u 



A 



for all u G N m , 1 < j ' < m, and 1 < k < n with u < (3—A k . 

(b): Let y G W be a nonnegative real solution to ( |^.^| )-( f^7l i and 
x = Ey, i.e. x k '■= Tli U y[k-,u\ for 1 < k < n. Then x belongs to the 
integer hull H of {x G M. n : Ax = b,x > 0}. Therefore a certificate 
of b ^ Ax for any x G N n is obtained as soon as any subsystem of 
(g^-(g7?]) has no solution y G W. 

Moreover, any nonnegative vector y G M p is a solution to ^^^ - ^fty 
if and only if y G V (with V as in h4-3\ )); and there is a one-to-one 
correspondance between the vertices of H and V. 
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Proof, (a) It is straightforward to see that f l4.4p -( TC7l) is just an explicit 
description of (13.1 7p in Theorem [71 The right- hand- side terms of (14.41) - 
(14.71) are the function b z —0 z evaluated at 2 6 N m with z < (3 and 
Ylij z j equal to 1, 2, 3, etc., until the maximal degree Yljftj- Each 
vector A^ z > G N p is then the transpose of the [zj-row of the product 
A© calculated in (I3.19p . so that 

A {z) [k,u} := (Ae)[z;(k,u)} = {{u+A k ) z -u z ) G N, 

for all appropriated indices z, k, and u. Evaluating (I3.17P at z = 
yields the trivial identity = Oy (recall Remark [5]), because we are 
using the convention that 0° = 1 and W = for every w^O. While 
evaluating (13.171) at z = ej, yields the linear constraints 

n 

bj = ^Aj. k y[ k > u ] V 1 < j < m, 

k=l ueN m , 
u</3-A k 

because w ej = Wj for all w G Z m . Constraints (14.41) are easily deduced 
after defining the vectors A^ G N p by A^ e ^[k,u] := Aj ]k for all u and 
1 < k < n. We also have that b = Ax with x = Ey and E e N" ,xp the 
matrix given in (14.31) . e.g. x k = ^2 u y[k,u] for 1 < k < n. 

(b) Let V E W be the polytope given in (j4\2j) . and H C R n be the 
integer hull of {x G R n : Ax = b, x > 0}. From Theorem H] and the 
proof of Theorem [7J any nonnegative vector y G MP is a solution to 
(I3.17P if and only if y G V, because 0y = b; and we have already 
stated in the proof of (a) that (I4.4l) - (l4.7p is completely equivalent to 
(I3.17p . We also deduced that x = Ey is a solution to Ax = b for 
any feasible solution y G V. Conversely, given any nonnegative integer 
solution x G N n to Ax = b and working as in the proof of Theorem HJ 
we can associate with x a nonnegative integer vector y(x) G N p such 
that 6y(x) = b and x = Ey(x); hence y(x) G V . 

Let {y^ G N p } be the vertices of V . Every y^ has nonnegative 
integer entries, because 6 is totally unimodular and the vector b in 
(14. II) has integer entries. Each vector Ey^ G W 1 has nonnegative 
integer entries and satisfies AEy^ = b, so that Ey^ G H. Any feasible 
solution y G V can be expressed as the sum y = Y^e&Y with = 
1 and ^£ > for each /. Hence, 

Ey = Y&Ey [e] with Ey [e] G H Ml, 

i 

and so x := Ey G H. Conversely, let {x 1 ^ G N n } be the vertices of 
H. We claim that every y(x^) is a vertex of V; otherwise if y(x^) = 
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£*&y [fl for 0<&<1, then 

X W = Ey(x [c] ) = J2& E y W with Ey w G if, 

is a combination of vertices (i£y™) of H, in contradiction with the fact 
that itself is a vertex of if. □ 

4.1. A hierarchy of LP-relaxations. Let \z\ : = 21+ • • ■ +z m for ev- 
ery z G N m . Consider the following integer and linear programming 
problems, for 1 < I < 

(4.8) J. := min { c'x : Ax = 6 , x G N n }; 

(4.9) J, := mm|c£y: ^ N m , 1 < |z| < i~ J" 
We easily have the following result. 

Corollary 10. Let (3, b G N m and A G N mxn 6e stzc/j tfiat (3>b and 
P > Ak for every index 1 < k < n. Let J, and Jp be as in fl^.<5| ) and 
h4-ty respectively. 

Then : J\ < J 2 < ■ ■ ■ < J\p\ and J\m = J,. 

Proof. It is obvious that Je < Je* whenever I < £*. Moreover : 
J. = min{c'x : x G H} and J\a\ = minjc'-Ey : y G V}, 

x y 

where H is the integer convex hull of {x G MJ 1 : Ax = b, x > 0} and V 
is given in ( 14.2R . Thus J, = c'x^ for some vertex x^ of if. Working as 
in the proof of Corollary [91 y(x^) is a vertex of V and 2^ = Ey(x^), 
so that 

J|^l < c'Ey(xW) = J.. 

On the other hand, Jigi = c'-Ey^ for some vertex y^ of V. Again, 
as in the proof of Corollary [9j Ey^ G ff , so that 

J. < c'EyM = J|^,. 

□ 

Hence the sequence of LP problems {Je}i<e<\p\ in (14.91) provides a 
monotone sequence of lower bounds for J. with finite convergence limit 
J\p\ — J»- Observe that all linear programs (14.91) have the same number 
p of variables (since y G MP) and an increasing number of constraints 
starting from m when £ = 1 to s = Ylj(0j J rl) when I = \/3\. The most 
interesting case happens when we take b = f3 > Ak (for every index 
1 < k < n). 
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Of course p (< ns) is large; but in principle the LP problems (14.91) are 
amenable to computation for reasonable values of I by using column 
generation techniques (in order to avoid handling the full vector y G 
W). However the vectors G N p may contain very large values 
when \z\ = i is large, and so some numerical ill-conditioned cases are 
likely to appear for large values of I. 

Therefore it is not completely clear whether the hierarchy of linear 
relaxations presented in Corollary [9] is useful from a computational 
point of view, but we think this hierarchy provides some new insights 
into integer programming problems from a theoretical point of view. 

4.2. A polynomial certificate. Let (3, b G N m and A G N mxn be 

such that (3 > b and (3 > A^ for every index 1 < k < n. Consider 
the matrices A G N sxs and 9 G Z sxp given in ([311- (D of Defini- 
tion [TJ Lemma [2] states that A is invertible, and so Theorem H] implies 
that existence of a nonnegative integer solution x G W 1 to Ax = b is 
equivalent to the existence of a nonnegative real solution y G MP to 

(4.10) AGy = Ab and y > 0, 

where b G V is given in (13 . 12[) or (14. ip . Moreover it was already stated 
in (I3.20p that the entries of the product Ab can be indexed with the 
nonnegative vector z G N m with z < (3, and that they have the simple 
representation : 

(4.11) (Ab)[z] = (^-0 2 ) = | ° = £ 

using the conventions 0° = 1 and W = for every w ^ 0. 

To obtain a polynomial certificate we apply Farkas lemma to the 
linear programming problem (14.101) . Consider the polyhedral cone 
Cfs C W defined by 

(4.12) Op := { £ G R s : (A6)'f > }. 

With every element £ = (£ 2 ) G Cb, we can associate a specific polyno- 
mial p^ G M[mi, . . . , u m ] defined by: 



(4.13) u^p ( (u):= Yl 



u 



zeN m , 

z^0,z</3 



and notice that in view of (14. lip . p^(b) = £'Ab. 

Theorem 11. Let f3, b G N m and A = [A l \ ■ ■ ■ \A n ] G N mxn be such that 
(3 > b and f3 > A k for every index 1 < k < n. Consider the polyhedral 
cone Cp C M s given in j^-lty and the following two statements : 



(a): The system Ax = b has an integer solution x G N r ' 
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(b): pz(b) < for some £ G Cp and p^ G M[u] as in M-13fy . . 
Then one and only one of statements (a) or (b) is true. 

Proof. The fact that A is invertible (because of Lemma [2]) and Theo- 
rem H] imply that the system Ax = b has a nonnegative integer solution 
x G N n if and only if (I4.10p holds for some real nonnegative vector 
y G MP. Farkas lemma implies that ( 14.101) holds if and only if £'Ab > 
for every £ G Cp] however from ( 14. lip . 

0<£'Ab = = 

2GN m , 
z^0,z<j3 



the desired result. 



□ 



It is interesting to compare Theorem [IT] with the standard Farkas 
lemma for (continuous) linear systems. The continuous Farkas lemma 
provides a linear certificate of the form uj'b < for some u G M. m that 
satisfies A'uj > 0; in contrast (the discrete) Theorem [11] provides a 
non-linear polynomial certificate p^(b) < 0. 

Example 12. Let A := [3,4] G N 2 , (3 = 5, and < b < 5; i.e. the 

Frobenius equation 3x\ + 4x 2 = b with Xi,x 2 G N. Then s = 6, 



A 



1111 i i 

1 2 3 4 5 

1 4 9 16 25 

1 8 27 64 125 

1 16 81 256 625 

1 32 243 1024 3125 



and 



-10 0-10 
0-10 0-1 

0-10 

1 
10 10 
10 1 



The pair ff^TW - ff^Tl ) reads as follows : 



- 











o ■ 


3 


3 


3 


4 


4 


9 


15 


21 


16 


24 


27 


63 


117 


64 


124 


81 


255 


609 


256 


624 


243 


1023 


3093 


1024 


3124 



y[i,o] 
i/[M] 

J/[l,2] 

y[2,o] 
y% i] 



o- 

6 
6 5 



Notice that A0 = [°] rai/i S 1 G N 5x5 an invertible matrix. Whence 
there is a nonnegative solution y > if and only if 

y (b) = S-%b 2 ,b 3 ,b\b B ]' > 0. 

We can easily verify that y(5) = [0, — 1, 0, 1, 1] ^ 0. We also have that 
y(b) >0forb = 0, 3, 4; and ^a* y(6) ^ /or 6 = 1, 2, 5. 

On the other hand, we can use the transpose A' instead of A in 
equations (I4.10p to (14.121) . Thus existence of a nonnegative integer 
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solution x G N n to Ax = b is equivalent to existence of a nonnegative 
real solution y G MP to 

A'0y = A'b, where (A'b)[z] = (z b - 1). 

Consider the polyhedral cone C% C M s defined by 

Q:={(GK S : 0'A£ > 0}, 

and with every £ = (£ z ) G CS, associate the exponential-like function 

(4.14) £ (*"-!)■ 

z€N m ,z</3 

Observe that ft(b) = £'A'b. By Farkas lemma, one and only one of the 
following statements holds : 

(a): The system Ax = z has an integer solution x G N n . 
(b*): f((b) < for some £ G (7g and / ? as in ( 14.14ft . 
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